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P Skizze Algorithmus 


1. Zufällige Initialisierung der Individuen 
2. Bewertung der Individuen 
3. Loop: Erzeugen neue Generation: 
(a) Selektion: Auswahl für nächsten Schritt 
(b) Anwendung genetische Operatoren auf 
selektierte Individuen 
(c) Bewertung der neuen Individuen 
(d) Erzeugen nächste Generation aus neuen und 
alten Individuen 
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P Genetische Algorithmen 


e Suche einer optimalen Lösung für ein Problem 

e Inspiriert von der Natur 

e Gene sind Sequenzen von Zahlen oder Symbolen 
(binär, Integer, Floating-Point) 

e Aus einem Gen wird eine potentielle Lösung für 
das Problem aufgebaut: Phänotyp 

e Algorithmus arbeitet mit Population einer 
bestimmten Größe 

e Fitness: Bewertungsfunktion wird optimiert 

e Maximierung oder Minimierung möglich 

¢ Individuum: Genotyp + Phanotyp 
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P Implementation in PGApack 


e PGApack: Parallel Genetic Algorithm Package 

e Open Source genetic algorithm implementation 

e Originally from David Levine, Argonne National 
Laboratory 

« Parallel execution via MPI (Message Passing 
Interface) 

e Construction kit for your experiments 

e I’m maintaining a fork since 2017 

e Current version implements differential evolution 
and other variants 

e PGApy is a python-wrapper (since 2005) 
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Nur fur Floating-Point Gene! [SP97, PSLO5] e Einfache „Box Constraints“ 

(a) Selektion: Nimm der Reihe nach ¢ Minimum und Maximum von jedem Parameter (=Gen) 
jedes Individuum k e z.B. bei Längen: Sind nicht negativ 

(b) Wähle zufällig drei Individuen, ¢ ... und meist ist Größenordnung bekannt 
alle verschieden und verschieden von k (mm, m, oder km) 

(c) Bilde Vektor-Differenz aus zwei dieser e Von praktisch allen genetischen Algorithmen 
Individuen und addiere Differenz zum dritten unterstützt 

(d) Optional: Crossover des Resultats mit k e Komplizierte Bedingungen: mehrere Zielfunktionen 

(e) nächste Generation: Das neue Individuum e Eval-Funktion hat mehrere Rückgabewerte 


ersetzt das alte wenn es besser ist 
Numerische Optimierung — auch für Antennen! 
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e Traditionell: Nur eine Zielfunktion 
(objective function) 
e Wird maximiert oder minimiert 
e Mehrere Zielfunktion: 
e für Randbedingungen (Constraints) 
e mehrere gleichzeitige Ziele 
— Randbedingungen erschweren das 
Optimierungsproblem 
— Ziele können sich widersprechen! 
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m: Beispiel 
e Vier Zahnräder 
e Vorgegebenes Übersetzungsverhältnis 
e Randbedingung: Zahnräder müssen sich gleich- 
mäßig abnutzen 
— Anzahl der Zähne von je einem kleinen und einem 
großen Rad ohne gemeinsamen Faktoren (GGT 1) 
e 17 < Anzahl Zähne < 90 
e Verhältnis 950:150 = 19:3 





Quelle: Wikimedia Commons, CC-by-sa 4.0 
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© sospir Annan 


e Lösung ohne Randbedingung: 17, 18, 38, 51 e Man kann alle Ziele in eine Funktion 

e Aber: 51 ist Vielfaches von 17 „vermantschgern” 

e und GGT (51, 18) = 3 und GGT (38, 18) = 2 e Das legt aber die Gewichtung von vorneherein fest 

e Lösung mit Randbedingung: 23, 49, 83, 86 « Vergleiche Vortrag letztes Jahr 

e Abweichung vom vorgegebenen Verhältnis: 0.00467% e Daher: Mehrere Ziele in einem Optimierungslauf 

e gefunden mit PGApy e Ziele werden sich typischerweise widersprechen 

e Differential Evolution, gerundet auf ganze Zahlen e Wir wollen also möglichst alle „guten“ Lösungen 
— examples/gears.py finden 

e Neues Project gearbox mit mehr mechanischen e Und später entscheiden welche Kompromisse wir 


Details eingehen 
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e Algorithmus von 2002 (!) Kalyanmoy Deb et. al [DPAMO2] 





e Non-dominated Sorting Genetic Algorithm 
e Sortiert nach Dominanz-Beziehung 
e Pareto-Front 
e Alle nicht-dominierten Lösungen entfernen 
e Zweite Front von nicht-dominierten Lösungen usw. 
« Letzte Front wird nur teilweise übernommen 
e Lösungen in Bereichen wo es wenige Lösungen 
gibt bevorzugt 
— Crowding Metric 
e Nach Konvergenz bleibt typisch eine Front 
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e Eine Lösung dominiert eine andere Lösung wenn 
e ... sie in mindestens einem Kriterium besser ... 
e... und in keinem anderen Kriterium schlechter ist 
« Alle nicht-dominierten Lösungen sind Pareto-Optimal 
e ... und bilden eine Pareto-Front (graphisch) 
e nach dem italienischen Wissenschaftler 
Vilfredo Pareto 
¢ ... der das Konzept der Pareto-Optimalität 
definiert hat 
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e Simulation mit NEC (numerical electromagnetics 
code) 

e Stammt aus den 80ern mit Lochkarten 

e Rewrite in C/C++ „nec2c“ 

e nec2c: Übersetzt ein NEC file in Tabellen-Ausgabe 

e xnec2: Nimmt ein NEC-File und zeigt Grafik 

e xnecview nimmt nec2c Ausgabe und zeigt Grafik 

e |.. alle in Debian gepackaged 


Antennen-Optimierung mit genetischen Algorithmen 
ist schon in den 90er Jahren publiziert [Lin97] aber 
vielleicht nicht genutzt weil bis 2015 patentiert? [AL95] 
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e Stehwellenverhaltnis VSWR 
(Voltage Standing Wave Ratio) 
e Wenn eine Antenne gut angepasst ist wird alles 


abgestrahlt 
e Wenn nicht wird ein Teil reflektiert 
— Stehwelle U U 
SWR = IV FR 
UVI-|UR| 


e Ist 1 wenn die Antenne gut angepasst ist 
e Gute Werte < 2 
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e Gewinn in Hauptrichtung 

e ... geteilt durch (max.) Gewinn 
in Ruckwartsrichtung 

e In dB 

e Intuitiv würde man annehmen das ist bei hohem 
Gewinn automatisch gegeben 

e... aber Antennen mit hohem Gewinn haben oft 
ausgeprägte „Keulen“ in der Gegenrichtung 
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e Gewinn gegenüber einem idealen Kugelstrahler 
e “Isotropic Radiator” 

e Kann man nicht bauen aber ist gut als Referenz 
e Einheit: dBi (Dezi-Bel über isotropic radiator) 

« Eine Antenne mit Richtwirkung hat Gewinn > 0dB 
e z.B. ein \/2-Dipol hat 2.15 dBi Gewinn 

e Gewinn gegenüber Dipol: dBd 

e Umrechnung dBi — dBd: subtrahiere 2.15 

° X dB entspricht einem Faktor 10*/1° 

e 3dB entspricht Faktor 2 

e 10dB entspricht Faktor 10 
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+ 4 Parameter 
e Coding als Floating-Point 
— Differential Evolution 
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Lange des zweiten Elements 
Abstand der Elemente 
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Eine Randbedingung: VSWR < 1.8 
Zwei Zielfunktionen: Gewinn und Vorwärts/Rückwärts 
Gen der Lange 4: 


e Lange des Falt-Dipols 
e Breite des Falt-Dipols 


Box-Constraints fur alle Parameter 





Two element antenna 


Forward/Backward ratio (dB) 





6.2 6.3 6.4 6.5 6.6 6.7 
Gain dBi 
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Two element, both fronts 


Forward/Backward ratio (dB) 
H 
` 


6.2 6.3 6.4 6.5 6.6 
Gain dBi 
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6.7 


6.8 
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e Publikation 2014 [DJ14] [JD14] 
e Funktioniert mit mehr Zielfunktionen 
Größe der Pareto-Front steigt exponentiell mit An- 
zahl Zielfunktionen! 
. macht NSGA-II ineffizient 
NSGA-III: Beschränkung auf wenige Punkte 
e Systematische Generierung von Ziel-Punkten 
e Nur die Crowding-Metrik wurde ersetzt! 
e Implementiert in PGAPack set 2021 
e Bis jetzt nur synthetische (Test-)Probleme gelöst 
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C3-DTLZ4 





e PGApack: 
github.com/schlatterbeck/pgapack 

e PGApy: 
github.com/schlatterbeck/pgapy 

e PyNEC: 
pypi.org/project/PyNEC/ 

« Antenna-Optimizer: 
github.com/schlatterbeck/antenna-optimizer 

e Gearbox: 
github.com/schlatterbeck/gearbox 
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